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Ce papier présente la notion de signal de 
transition de phase pour un système de taille 
finie. Il se concentre sur le rôle des contraintes 
physiques sur ces signaux et la robustesse quant 
au changement d'ensemble statistique. Des 
résultats obtenus avec un modèle de gaz sur 
réseau sont présentés en analogie avec les 
observations effectuées en collisions d'ions 
lourds aux énergies intermédiaires. 



1 Introduction 

Les transitions de phase dans les systèmes finis ont été 
étudiées dans de nombreux cadres théoriques. Des ap- 
proches simples comme le modèle de gaz sur réseau per- 
mettent d'étudier certaines propriétés comme la coexis- 
tence de phases ou la criticité et décrire de façon qual- 
itative les processus physiques comme la multifragmen- 
tation des noyaux atomiques ou des agrégats d'atomes. 
La description de la transition de phase dans ces sys- 
tèmes diffère de la phénoménologie classique que nous 
connaissons bien. Pour ces systèmes de taille finie, les 
fluctuations des observables sont importantes, les fonc- 
tions de partition contenant un nombre fini de termes 
ne permettent pas aux potentiels thermodynamiques de 
diverger, le rôle des interfaces n'est plus négligeable et 
surtout les ensembles statistiques ne sont pas équiva- 
lents. Une nouvelle description des transitions de phase 
pour ces sytèmes doit être formulée. Après une présen- 
tation du modèle de gaz sur réseau et quelques rappels 
de physique statistique, je présenterai les effets des con- 
traintes physiques appliquées au système, c'est-a-dire le 
choix d'un ensemble statistique, sur deux signaux de 
transition de phase. 



2 Le modèle Lattice-gas 

Le modèle du gaz sur réseau est une extension du modèle 
d'Ising. Ce modèle suppose un réseau multidimensionnel 
(3 dimensions dans notre cas) comprenant N sites carac- 
térisés par un nombre d'occupation rij (absence (0) ou 
présence (1) d'une particule sur le site i) ainsi qu'une 
impulsion pl. La donnée des N nombres d'occupation 
et impulsions associées définit une configuration micro- 
scopique {n}. Le Hamiltonien que nous utilisons dans 



nos simulations s'écrit de la façon suivante 
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H({n}) 



fc=l 



La première somme traduit une interaction attractive 
constante entre plus proches voisins sur le réseau (6 au 
maximum en 3 dimensions) . Le paramètre e peut être vu 
comme la profondeur du puits de potentiel d'une inter- 
action générique à courte portée. Cette première somme 
est responsable de la cohésion du système et joue un rôle 
de première importance. La deuxième somme prend en 
compte l'agitation du système à travers l'impulsion des 
particules du réseau donnant ainsi des degrés de lib- 
erté supplémentaires au système. Cette agitation con- 
tribue à la dissociation du système en fragments. Les 
fragments sont identifiés à l'intérieur de chaque config- 
uration microscopique par un algorithme de reconnais- 
sance du type Coniglio-Klein fJJ. L'observation du sys- 
tème consiste ensuite à évaluer des observables perti- 
nentes calculées sur un état microscopique du système. 



3 Ensembles statistiques 

Les contraintes physiques extérieures sur le système 
étudié sont appliquées par le biais de multiplicateurs 
de Lagrange. L'observable Ai (variable extensive) asso- 
ciée à A; (paramètre intensif) se retrouve alors distribuée 
dans un intervalle. Par exemple le multiplicateur de La- 
grange associé au volume est la pression, celui associé 
au nombre de particules le potentiel chimique ou en- 
core celui associé à l'énergie la température. Pour toute 
observable nous pouvons définir un paramètre de La- 
grange (conjugué canonique) qui permettra d'obtenir la 
distribution de cette observable à l'équilibre de Gibbs, 
c'est à dire celle maximisant l'entropie statistique sous 
contraintes. Cette méthode permet de traiter de façon 
équiprobable l'information dont on ne dispose pas, le 
seul biais introduit étant précisément l'information per- 
tinente. 

Un algorithme Monte-Carlo de type Metropolis [Sjj per- 
met d'échantilloner différents ensembles statistiques afin 
d'obtenir les distributions statistiques de microétats sat- 
isfaisant les contraintes. Il est possible d'ajouter à ces 
contraintes des lois de conservations (distribution de 
Dirac de l'observable conservée). Nous considérerons 
les ensembles grand-canonique dans lequel le potentiel 
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Lois de conservation 
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Distribution de probabilité 
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Tab. 1 - Tableau récapitulant les propriétés des ensembles statistiques utilisés. 



chimique et la température sont fixés (ie. nombre de par- 
ticules et énergie définis en valeur moyenne) ainsi que 
l'ensemble canonique où la masse est fixée et l'énergie 
est définie en valeur moyenne. Les distributions de prob- 
abilité pour chaque ensemble statistique sont rappelées 
dans le tableau 



4 Bimodalité 

Une fois les distributions statistiques de microétats 
obtenues, les projections de la densité d'état sur des 
directions d'observation jugées pertinentes permet d'é- 
tudier la physique associée au système. 

F. Gulminelli et P. Chomaz ont proposé |2J d'utiliser la 
topologie de la distribution des événements dans l'espace 
d'observation pour signer la transition de phase pour 
les systèmes de taille finie. Dans la zone de coexistence 
de phase une bimodalité (deux maxima) dans la distri- 
bution de probabilité d'une observable permet d'identi- 
fier sans ambigûité chacune des phases ; cette observable 
joue alors le rôle d'un paramètres d'ordre. La projection 
de la densité de probabilité dans des espaces d'obser- 
vation pertinents permet d'identifier l'appartenance des 
événements à chacune de ces phases. 

La densité de matière est le paramètre d'ordre naturel 
de la transition de phase liquide-gaz à la limite thermo- 
dynamique. Il est toutefois possible de trouver d'autres 
paramètres d'ordre corrélés à la densité ; en particulier 
on peut s'attendre à ce que la taille du plus gros des 
fragments {Abig) puisse servir de paramètre d'ordre. 

La figure Q représente la distribution des événements 
grand-canoniques (à la valeur critique du potentiel chim- 
ique /i = 3e) dans l'espace "taille totale du système (ie. 
densité) versus taille du plus gros fragment" et montre 
la corrélation entre ces deux observables. Cette figure 
montre que même en système de toute petite taille la 
coéxistence de phase est reconnaissable par la bimodal- 
ité de la distribution du paramètre d'ordre. La claire 
corrélation entre Au g et et la masse totale du système 
A to t démontre en plus que A^g peut être utilisé comme 
paramètre d'ordre. Cette observable est facile à manip- 
uler dans les simulations numériques et accessible assez 
facilement expérimentalement. La densité étant beau- 
coup plus difficile (voire impossible du fait de la diffi- 
culté d'accéder au volume) à mesurer, l'étude de cette 
observable est d'un grand intérêt pour faire des analogies 
ou comparer les résultats des modèles avec l'expérience. 
A basse température la distribution de ces observables 
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Biggest fragment size 

FiG. 1 - Densité de probabilité grand-canonique pro- 
jetée dans l'espace "nombre total de fragments" (ordon- 
nées) versus "taille du plus gros fragment" (abscisse) 
pour différentes températures (dans un réseau 6x6x6) 



est bimodale (la projection sur chaque axe révèle deux 
maxima). Les événements denses (grande masse totale) 
et ayant un gros plus gros fragment Au g sont identifiés 
comme appartenant à la phase de type liquide tandis 
que le deuxième lot d'événements (faible densité et pe- 
tite taille du plus gros fragment) correspond à la phase 
gazeuse. En élevant progressivement la température du 
système, ces deux lots d'événements se rapprochent l'un 
de l'autre pour fusionner à la température critique. En- 
fin pour les températures surcritiques la distribution de- 
vient monomodale (un maximum), une seule phase flu- 
ide subsiste, avec un plus gros fragment de petite taille. 
Ainsi le passage d'une distribution bimodale à une dis- 
tribution monomodale signe le passage par un point cri- 
tique correspondant à la transition de phase du second 
ordre. 

La figure (01 (graphique supérieure) schématise la bifur- 
cation caractéristique d'une transition de phase. Cette 
figure montre la distribution grand-canonique du plus 
gros fragment pour différentes températures. Le pas- 
sage de la distribution bimodale de ce paramètre d'or- 
dre à une distribution monomodale est bien localisée et 
permet d'identifier le point critique. Le graphique in- 
férieur sur cette même figure présente les distributions 
du plus gros fragment dans la même gamme de tempéra- 
tures pour l'ensemble canonique. Dans cet ensemble la 
masse a été fixée, c'est-à-dire le nombre de sites occupés 
est maintenu constant, alors qu'il était libre de fluctuer 
(sous la contrainte du potentiel chimique) dans l'ensem- 
ble grand-canonique. Si l'existence d'une loi de conserva- 
tion empêche évidement l'utilisation de la masse totale 
comme paramètre d'ordre, la bimodalité de la distribu- 
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Grand Canonical Transition In space Zbiq vs T 




Canonical Transition in spac& Zbig vs T 




FlG. 2 - Distribution de taille du plus gros fragment 
pour différentes températures dans l'ensemble grand- 
canonique (haut) et dans l'ensemble canonique (bas). 



tion de taille du plus gros fragment a aussi disparu dans 
cet ensemble. La contrainte sur l'observable "masse to- 
tale du système" a fondamentalement changé la topolo- 
gie de la distribution de A big à cause de la forte corréla- 
tion entre A big et A tot . Le plus gros fragment n'est plus 
un paramètre d'ordre pour cet ensemble et ne permet 
plus de caractériser la coéxistence de phase. 



Invariance 
d'échèlle 



par changement 



Il a été proposé de nombreux autres moyens de mettre en 
évidence une transition de phase. R. Botet et M. Plosza- 
jczak t 4| ont proposé par exemple d'utiliser la théorie des 
fluctuations universelles ou A — scaling. Cette théorie 
démontre que dans un phénomène critique la distribu- 
tion du paramètre d'ordre est invariante par changement 
d'échèlle du système. Cette théorie a des applications en 
physique nucléaire pour ce qui concerne la phénoménolo- 
gie de la multifragmentation des noyaux excités à des én- 
ergies comparables à leur énergie de liaison, si la taille 
du plus gros fragment est paramètre d'ordre de la tran- 



sition de phase associée à la multifragmentation et si 
celle-ci est interprétable comme un phénomène critique 
(ie. a lieu en proximité du point critique du diagramme 
de phase) . On attend dans ce cas que la forme de la dis- 
tribution du plus gros fragment (moyenant une trans- 
formation d'échèlle) reste la même lorsque la taille du 
système est modifiée traduisant le caractère universel 
de cette distribution. Le passage par un point critique 
se traduit par la vérification du A — scaling de la dis- 
tribution normalisée -P< m >(m) du paramètre d'ordre m 
pour différentes tailles de systèmes : 



< m > A P <m> (m) = $(za), z a = 



< m > A 



(1) 



où m* et < m > représentent la valeur la plus probable 
et la valeur moyenne du paramètre d'ordre et <& la fonc- 
tion d'échèlle. Dans une phase désordonnée les distribu- 
tions des paramètres d'ordre doivent vérifier A = 1, le 
passage par un point critique est signé par A = \. J.M. 
Carmona et al. ont montré que cette théorie est véri- 
fiée dans le modèle de gaz sur réseau. Cette technique 
est cependant difficle à mettre en oeuvre car la qualité 
de la loi d'échèlle (QJ dépend évidemment de la valeur 
de A qui peut à priori varier entre et 1 et n'est pas 
connue avant l'analyse. 

En ce qui concerne les données expérimentales, le nom- 
bre limité d'espèces nucléaires rend impossible une ap- 
plication directe de la méthode. Les propriétés d'invari- 
ance d'échèlle ont alors été étudiées pour une même colli- 
sion d'ions lourds en variant la violence de la collision 
et un changement de régime a été observé de A = 1 
vers A = \. L'ensemble statistique (tri expérimental) 
dans lequel est effectué cette analyse s'apparente à un 
ensemble gaussien intermédiaire entre un ensemble mi- 
crocanonique et un ensemble canonique. 

En notant P(A b i g ) la distribution du plus gros fragment, 
< Au g > et (TAug la moyenne statistique et la variance 
de cette distribution, l'équation peut être écrite de 
façon équivalente par l'ensemble des deux relations 



a Au P(A Ug ) = V(z P ms) = *( 
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= Constante 



(2) 



(3) 



L'équation (J2J exprime la transformation qui permet 
d'obtenir les distributions centrées ^{zpms) de la vari- 
able réduite zpms- Ceci constitue une première étape 
qui ne suppose pas la valeur de A. L'équation Q 
exprime une relation affine entre les logarithmes des 
valeurs moyennes et des variances des distributions de 
plus gros fragments 

ln(A big ) = Cte - Aln(a abig ) 

qui permet, si la loi d'échèlle l|2} est vérifiée, d'établir la 
valeur de A. 

Nous avons appliqué cette méthode d'analyse aux don- 
nées de notre modèle pour mieux comprendre ce signal. 
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FiG. 3 - Logarithme de la variance de la distribution de 
taille du plus gros fragment en fonction du logarithme 
de la valeur moyenne (figure de gauche) et fonctions 
de scaling (cf. texte) de ces distributions (droite) dans 
l'ensemble grand-canonique (haut) et canonique (bas). 

Similairement à la méthode expérimentale, nous avons 
gardé la taille du réseau constante (6x6x6, ce qui corre- 
spond à 108 particules dans l'ensemble canonique). 

La figure (01 montre les résultats obtenus dans les en- 
sembles grand-canonique (haut) et canonique (bas). La 
partie gauche de chaque sous-figure montre la relation 
entre le logarithme de la valeur moyenne (Zn(/tf)) et 
le logarithme de la variance (Zn(«2)) des distributions 
du plus gros fragment, chaque point correspond à une 
température (basses températures à droite, hautes à 
gauche), la partie droite de ces figures superpose les 
fonctions d'échelle #(zpms)- 

Dans l'ensemble grand-canonique, comme nous l'avons 
vu au paragraphe précédent, les distributions du plus 
gros fragment dans la zone de coexistence de phase sont 
bimodales. La variance globale a donc une signification 
qui n'est plus compatible avec la théorie des fluctuations 
universelles. Dans cet ensemble les fonctions d'échelle 
de ces distributions ne se superposent pas, il n'y a pas 
d'invariance par changement d'échèlle. La relation en- 
tre variance et moyenne est linéaire, aucun accident ne 
permet dans ce diagrame de déceler la localisation du 
point critique. Il est plus judicieux par exemple d'utiliser 
comme nous l'avons vu la disparition de la bimodalité 
pour localiser le point critique. Seules les distributions 
correspondant à des températures surcritiques semblent 
approximativement suivre une loi d'échèlle avec A sa 1. 

Dans l'ensemble canonique, les conclusions sont très 
différentes, la valeur moyenne de la taille du plus 
gros fragment diminue avec la température, la vari- 
ance augmente, passe par un maximum puis rediminue 
traduisant l'élargissement de la distribution de plus gros 
fragments observée au voisinage de la température cri- 



tique. Les fonctions d'échèlle sont regroupées en deux 
familles distinctes possédant une asymétrie propre. Les 
distributions correspondant aux basses températures 
sont inclinées vers la droite tandis que celles correspon- 
dant aux grandes températures sont penchées vers la 
gauche, au passage de la température critique cette dis- 
tribution du paramètre d'ordre s'aplatit avec une forme 
très large (fluctuation maximale). Assez loin de la tem- 
pérature critique, les distributions de chaque famille se 
ressemblent beaucoup mais ne suivent pas rigoureuse- 
ment une loi d'échèlle. Le A — scaling n'est pas observé 
bien que le point critique thermodynamique soit exploré 
par la simulation (point d'intersection des deux droites 
dans le diagramme de gauche). En effet, même si on 
accepte d'interpréter la ressemblence des distributions 
comme une loi d'échèlle déformée par les effets de taille 
finie, il est clair sur la figure Q que la relation entre 
In(Abig) et ln(<JA big ) n'est pas linéaire sauf peut-être 
pour les températures fortement surcritiques. Toutefois 
les interpolations linéaires de la figure © montrent que 
si l'on dispose d'un domaine de variation de tempéra- 
ture limité, la réduction de la variance due à la loi de 
conservation sur la masse totale peut être faussement 
interprétée comme un changement de régime de A = 1 
àA = |. 

6 Conclusion 

Dans cette contribution nous avons montré que les 
distributions des observables et les moments associés 
(moyenne, variance ...) dépendent de façon critique de 
l'ensemble statistique considéré (ie. les contraintes ap- 
pliquées à notre système). Ceci est spécialement impor- 
tant dans la recherche d'un paramètre d'ordre pour une 
transition de phase. Une observable peut être paramètre 
d'ordre pour un ensemble statistique et ne plus l'être 
pour une autre. En particulier notre étude montre qu'il 
est difficile de penser que Au g puisse être utilisé comme 
paramètre d'ordre dans les expériences de multifragmen- 
tation nucléaire étant donné sa corrélation avec la masse 
totale du système qui est fortement contrainte dans les 
expériences. Par contre des observables se basant sur 
l'asymétrie des distributions du plus gros fragment sem- 
blent très prometteuses. 
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